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Abstract. We present a simple way to derive the results of Diaconis and Fulman 
I ^arXiv: 1102.515 9' in terms of noncommutative symmetric functions. 

D 

^ ■ 1. Introduction 

In 121 13] , Diaconis and Fulman investigate a remarkable family of matrices intro- 
O ' duced by Holte ^ in his analysis of the process of "carries" in the addition of random 

integers in base b. 

^ ■ The aim of this note is to show that the results of [3] can be derived in a simple 

', and natural way within the formalism of noncommutative symmetric functions [8]. 

This is possible thanks to the following equivalent characterization of the "amazing 
matrix" P (Theorem 2.1 of |3j): 

The number of descents in successive 6-shufHes of n cards form a 
Markov chain on {0, 1, . . . , n — 1} with transition matrix P{i, j). 

Such random processes involving descents of permutations can usually be inter- 
preted in the descent algebra of the symmetric group. Here, it is only the number of 
descents which is involved, so that one can in fact compute in the Eulerian subalgebra. 



> 
OO 



0^ 

O I We assume that the reader is familiar with the notations of 

2. The Eulerian algebra 



X 



This is a commutative subalgebra of dimension n of the group algebra of the 



' symmetric group ©„, and in fact of its descent algebra S„. It was apparently first 



introduced in pLj under the name algebra of permutor^. It is spanned by the Eulerian 
idempotents, or, as well, by the sums of permutations having the same number of 
descents. 

It is easier to work with all symmetric groups at the same time, with the help of 
generating functions. Recall that the algebra of noncommutative symmetric functions 
Sym is endowed with an internal product *, for which each homogeneous component 
Sym„ is anti-isomorphic to E„. 

Recall also the following definitions from [8]. We denote by at or (7t{A) the gener- 
ating series of the complete symmetric functions Sn'- 

(1) at{A) = Y,t''S^{A). 

n>0 



^ A self-contained and elementary presentation of the main results of can be found in [5] . 
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\k] 

The Eulerian idempotents En are the homogenous components of degree n in the 
series E^'^^ defined by 

(2) = ^x'^^WlA). 

fc>0 

We have 

n 

(3) * ijH = , and ^i?W = 5„, 

k=l 

SO that the En^ span a commutative n-dimensional *-subalgebra of Sym„, denoted 
by Sn and called the Eulerian subalgebra. 

The noncommutative Eulerian polynomials are defined by 

n n 

(4) An{t) = E ( E ^0 = ^ A{n,k)t', 

k=l \I\=n k=l 

e(i)=k 

where Ri is the ribbon basis. The following facts can be found (up to a few misprints) 
in [S]. The generating series of the An{t) is 

(5) A{t) := E Mt) = (1 - t) (1 - t ^1-*)-' . 

n>0 

Let Al{t) = (1 - t)-" Then, 

/ \ ^(J) 

(6) ^= E = E T37 

n>0 7 ^ 

This last formula can also be written in the form 

(7) A*{t)= E (yZ-iX (Si + S, + Ss + ---)' 

k>0 ^ ^ 

or 



Let iSt*^' = ai{A)'' be the coefficient of t'^ in this series. In degree n, 



I^nl{I)<k ^ ^ i=l 



(9) 

This is another basis of £n- Expanding the factors (1 — in the right-hand side 
of ([H]) by the binomial theorem, and taking the coefficient of in the term of weight 
n in both sides, we get 

i 



(10) E rtMA("'^-^)- 



i=0 
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Conversely, 



fc>0 



SO that 

(12) A{n,p)= E (-ir t ^) • 

i=o ^ ^ 

The expansion of the En'^ on the basis A(n, i), which is a noncommutative analog of 
Worpitzky's identity (see [7] or [10]) is 

(13) 5: (^ + "-MA(n,.). 

fc=l i=\ ^ ' 

Indeed, when x is a positive integer A^, 

n 

(14) = ^«(^^) = E 

k=l H=n 

where Fj are the fundamental quasi-symmetric functions, and for a composition / = 
(zi, . . . , v) of n. 



(15) Fj{N) 



N + n — r 
n 



3. The ^-shuffle process 

For a positive integer 6, the 6-shuffle permutations in (3„ are the inverses of the 
permutations with at most h—1 descents. Thus, the 6-shufHe operator can be identi- 
fied with {i.e., with ^-multiplication by S^). It belongs to the Eulerian algebra, 
so that it preserves it, and it makes sense to compute its matrix in the basis A(n, k). 
Note that since £n is commutative, it does not matter whether we multiply on the 
right or on the left. 

Note that the 6-shuffle process is an example of what Stanley has called the QS- 
distribution ^JJ. It is the probability distribution on permutations derived by as- 
signing probability to the first h positive integers, see [S| for a simplified version. 

Summarizing, we want to compute the coefficients Pij{b) defined by 

n 

(16) S^UA{n,J) = Y,P^Ab)Mn,^). 

i=l 

From ([9]), it is clear that 

(17) * = ^[f^] 
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SO that, using (|T2l) . we obtain 



^[fl * A{n, j) = ^(-l)" + ^) SMj - r)] 

.1 — n \ / z — n V / 



:i8) '"=° 



r=0 ^ ^ A:=0 

The coefficient of A(n, z) in this expression is therefore 

(19) ^..w-B-irC:')^^"^''"') 

These are the coefficients of the amazing matrix (up to a shift of 1 on the indices 
and a global normalization factor fc" so as the probabilities sum up to 1). 

Since the En^ form a basis of orthogonal idempotents in Sn, it is reasonable to 
introduce a scalar product such that 

(20) (4^1|4^'1) = 5.,- 

Then, the 6-shuffie operator is self-adjoint. Its orthonormal basis of eigenvectors is 
clearly En^ (with eigenvalues b''). 

In terms of coordinates, since we are working in the non-orthogonal basis A{n,i), 
its right eigenvector of eigenvalue is the column vector whose ith component is the 
coefficient of eH^ on A{n,i), that is, the coefficient of in thanks to ( fT3l) . 

By duality, its left eigenvector associated with the eigenvalue 6* is the row vector 
whose jth component is 

(21) {A{n,j)\Ell) = j^i-^y ^ ^) ij - 

This is precisely the Foulkes character table (up to indexation, the Frobenius char- 
acteristic of x"''^ is the commutative imag^ of A(n, n — k)). 

4. Other examples 

4.1. Determinant of the Foulkes character table. This is the determinant of 
the matrix F 

(22) F{i,j) = {Ain,i),El^^) z,j = l,...,n. 
Because of the triangularity property 



i—r\ 
n ' 



(23) Mn,^) = Sli^ + j2i-^y(''^^)sl 

r=l V ^ / 

we have as well 

(24) det F = det G where G{i, j) = {SH^ , ^jf 1) = 



These commutative symmetric functions have been studied in [5]. 
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a Vandermonde determinant which evaluates to n!(n — 1)! ■ ■ ■ 2!l!. 
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4.2. Descents of fe^'-riffle shuffles. Recall from [5] that FQSym is an algebra based 
on all permutations and that it has two bases 

(25) Yl «' = F--^ 

std(to)=(T 

which are mutually adjoint for its natural scalar product 

(26) {F^,Gr) = 6„^r. 

Under the embedding of Sym into FQSym, the ^''-shuffle operator is 

(27) {S^Y = ST^= Yl F.. 

(T fe*"— shuffle 

The generating function of ^''-shuffle by number of descents is therefore its scalar 
product in FQSym with the noncommutative Eulerian polynomial 



(28) Anit) = Y ^'A(n, k)= J2 f^'^^^Gr. 

k=l tG6„ 

Recall that 

n 

(29) An{t) = {l-tr+'Y.^'Sl^^ 

k=l 

so that 

n 

(30) {sr,Mt)) = ii-tr^'J2t'isr,si^^) 

k=l 

Now, when one factor P of a scalar product {P, Q) in FQSym is in Sym, one has 
(P, Q) = {p, Q) where p = P_ is the commutative image of P in QSym, and the 
bracket is now the duality between Sym and QSym. Furthermore, when p in in 
Sym, then, the scalar product reduces to {p, q) , where q = Q is the commutative 
image of Q in Sym, and the bracket is now the ordinary scalar product of symmetric 
functions (see [5j). Thus, 



(31) (5f Ui") = {K{b'X),KikX)) = K{u-k) 

(A-ring notation) and we are done: 

n . 

(32) {€\Mt)) = {i-tT^'Y.^H 

k=l ^ 



b^'k + n-l 
n 



Ifk + n-l 
n 
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